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Perturbations of the linearized vacuum Einstein equations on a null cone in the Bondi-Sachs 
formulation of General Relativity can be derived from a single master function with spin weight 
two, which is related to the Weyl scalar *&o, and which is determined by a simple wave equation. 
Utilizing a standard spin representation of the tensors on a sphere and two different approaches 
to solve the master equation, we are able to determine two simple and explicitly time-dependent 
solutions. Both solutions, of which one is asymptotically flat, comply with the regularity conditions 
at the vertex of the null cone. For the asymptotically flat solution we calculate the corresponding 
linearized perturbations, describing all multipoles of spin-2 waves that propagate on a Minkowskian 
background spacetime. We also analyze the asymptotic behavior of this solution at null infinity 
using a Penrose compactification, and calculate the Weyl scalar, ^4. Because of its simplicity, the 
asymptotically flat solution presented here is ideally suited for testbed calculations in the Bondi- 
Sachs formulation of numerical relativity. It may be considered as a sibling of the well-known 
Teukolsky-Rinne solutions, on spacelike hyp ersurf aces, for a metric adapted to null hypersurfaces. 

PACS numbers: 04.20.-q, 04.20.Jb, 04.25.-g, 04.25.D- 
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I. INTRODUCTION 

Exact solutions of the Einstein field equations provide 
a deeper insight into the classical theory of General Rela- 
tivity. In this article, we present an exact time-dependent 
global solution on a null cone describing all multipoles 
of linearized spin-2 fields propagating on a Minkowskian 
background space by using the Bondi-Sachs formulation 
of General Relativity. The solution is given as a 
spectral series with respect to spin-2 harmonics, where 
the coefficients are simple rational expressions of the time 
and radial coordinate. Therefore, we refer to it as SPIN- 
2. It is an ideal textbook solution allowing one to demon- 
strate, when working in the Bondi-Sachs frame work of 
General Relativity, important features, such as the reg- 
ularity conditions at the vertices, the commonly used 
3— formalism, and the subtleties at null infinity. In ad- 
dition, since it describes all radiation multipoles, SPIN-2 
is also well suited as a testbed solution for numerical rel- 
ativity, when the Einstein field equations are solved in 
a Bondi-Sachs framework, (see [4J for a review). Thus 
it might be considered as a sibling, in null coordinates, 
of the well-known Teukolsky-Rinne solutions p], [|| em- 
ployed in numerical relativity using the (3 + 1) formula- 
tion of General Relativity. 

Despite its simplicity, this is the first time that a regu- 
lar and asymptotically flat solution of the linearized vac- 
uum Einstein equation has been reported for all multi- 
poles in Bondi-Sachs formulation. Linearized solutions 
were first discussed qualitatively by Bondi et al. [l[ , who 
gave an asymptotic vacuum solution in terms of inverse 
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powers of an areal distance coordinate r of an axisymmet- 
ric metric with a hypersurface orthogonal Killing vector. 
As their solution was given by coefficients of a series of 
r~ n , (n > 0), it is not regular at the vertices of the null 
cones. Winicour proposed a Newtonian approach to 
the Bondi-Sachs formulation of General Relativity which 
is related to a post-Minkowskian expansion of the Bondi- 
Sachs metric. In particular, he pointed out the necessity 
of imposing regular boundary conditions at the vertex of 
a freely falling Fermi observer in the background space- 
time, and he introduced spin-0 potentials (l0| to solve 
for the perturbations. Axisymmetric linearized solutions 
were revisited by Papadopoulos [III EH > wn0 use d Wini- 
cour's idea of the spin-0 potentials to find solutions of 
linearized vacuum perturbations. Papadopoulos's algo- 
rithm was later generalized by Lehner jl3l [l4| to three 
dimensions. 

The spin-0 potential approach to find a radiative 
I— multipolc can be summarized in three basic steps: 
First; one guesses a regular solution at the vertex for a 
monopole scalar field that obeys the flat space wave equa- 
tion. Second; one applies 71— times the z— translation op- 
erator expressed in outgoing polar null coordinates [39| to 
the monopole solution to find a n— multipolc of the scalar 
wave equation. These multipoles are also a solution of the 
scalar wave equation, because the z— translation operator 
commutes with the axially symmetric angular momen- 
tum operator. Finally; one finds the Bondi-Sachs metric 
functions by applying the 9— operator [l(], EH to the 
?i— multipoles solution of the scalar wave equation and 
integrates it to obtain the Bondi-Sachs metric functions. 
This elegant method has the disadvantage that it requires 
an infinite application of the z— translation operator for 
generating a solution for all n— multipoles. 

Linearized and quadratic perturbations, with respect 
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to a Minkowski background, were considered by [l6| . who 
presented their solutions in terms of Newman-Penrosc 
quantities [13] and, like Bondi et al. [l|, they gave only 
leading order terms of a r -1 expansion of these quanti- 
ties in the asymptotic regime. Bishop and collaborators 
[3 EH proposed a procedure to find linearized perturba- 
tions on a Schwarzschild background. However, in their 
approach "...some of the expressions get very complicated 
and in order to simplify the presentation we now special- 
ize to the case I — 2..." [l8| to obtain the solution is 
overly complicated, although the calculation of linearized 
perturbations is achievable by a simpler method. Never- 
theless, Reisswig et al. [2(j determined a I = 3 multipole 
using Bishop's approach. 

We solve the vacuum Einstein equations for the zeroth 
and first order terms of an expansion of the metric in 
terms of a measure of the deviation from spherical sym- 
metry. Thereby we assume that the null cones emanate 
from a Fermi observer [2l| following the timelike geodesic 
of the background spacetime. As the boundary condi- 
tions at the vertices are given by the regularity conditions 
of the Fermi observer [221 ] , we determine these boundary 
conditions for the general three dimensional case using 
spin weighted harmonics. Utilizing these boundary con- 
ditions, we then integrate the Einstein equations from 
the vertex to infinity employing two different approaches: 
one in which an asymptotically flat solution is obtained 
whereas in the other one the solution diverges exponen- 
tially. Applying the Penrose compactification [23| to the 
asymptotically flat solution, the Weyl scalar ^4 at null 
infinity is calculated with the formalism of 24 1 . 

This article is organized as follows: In Sec. HH we intro- 
duce the quasi spherical approximation of a Bondi-Sachs 
metric. In Sec. Mil we derive the necessary equations and 
framework to find solutions of the zeroth and first-order 
quasi-spherical approximation. In particular, we intro- 
duce a function, which we call master function, which 
is related with the Weyl scalar ^>q, and whose solution 
allows us to determine the linearized perturbation of a 
vacuum spacetime. We also present equations obtained 
with the 3— formalism that will be used to find a solution 
of the perturbations. In Sec. HVi we determine the gen- 
eral boundary conditions of the master function and the 
resulting boundary conditions for the perturbations in 
three dimensions using spin-spherical harmonics. In Sec. 
fVl two solutions of the master function are presented, 
one of which diverges exponentially as r tends to infinity 
whereas the other one is asymptotically finite. For the 
asymptotically finite solution of the master function, we 
calculate the resulting perturbations and the Weyl scalar 
p4 at null infinity in Sec. I VII adopting the formalism of 
[23 ]. Finally, our results are summarized in Sec. I VIII 

We use geometrized units in our calculations and the 
conventions of [HI for curvature quantities. For tensor 
indices, we use the Einstein sum convention, where small 
Latin indices (a, b, c, ...) take values 0,..., 3 and capital 
Roman indices (A, B, C, ...) have values 2, 3 correspond- 
ing to angles 9 and 0, respectively. Bondi-Sachs coordi- 



nates are denoted by x a , conformal Bondi-Sachs coordi- 
nates by x a , and conformal incrtial coordinates by x a . 
A quantity calculated in conformal Bondi-Sachs coordi- 
nates is denoted by a hat over the respective quantity, one 
in conformal coordinates by a tilde, and a complex con- 
jugated quantity by an over-bar. The expression A=B 
means that A — B = 0(e 2 ) in a quasi-spherical approxi- 
mation. 



II. QUASI-SPHERICAL APPROXIMATION OF 
THE BONDI-SACHS METRIC 

Consider a smooth one-parameter family of metrics 
g a b{^) at null cones in a vacuum spacetime, where e is a 
parameter that measures deviations from spherical sym- 
metry, and when e = we recover the metric of a spher- 
ically symmetric spacetime. This spherically symmet- 
ric spacetime is referred to as background spacetime and 
contains a unique geodesic, the central geodesic, which 
traces the centers of symmetry. The metrics g a b(£) are 
expressed in terms of Bondi-Sachs coordinates, x a , [lH3- 
Given the central timelike geodesic of the background 
spacetime, the coordinate x° = u is the proper time 
along the geodesic and is constant along outgoing null 
cones for points that are not on this geodesic. The two 
coordinates x A = (x 2 , x 3 ) are angular coordinates pa- 
rameterizing spheres centered on points on the timclikc 
geodesic. Finally, the coordinate x 1 — r is an areal dis- 
tance coordinate, such that surfaces dr = = du have 
the area 4irr 2 . The line element of the metric is given by 

ds 2 = _ e »( s )-H/»( S ) du 2 _ 2e 2 ^dudr 

+r 2 h AB (e)[dx A - U A (e)du][dx B - U B (e)du] (1) 

where h AC = h C B = S A B and |/mb|,t- = = |/iab|,« 
Q . The latter condition on the metric 2-tensor Hab is a 
consequence of the requirement of the radial coordinate 
to be an areal distance coordinate. This implies that Kab 
has only two independent degrees of freedom to describe 
the geometry of the 2-surfaces dr = = du. The metric 
functions in ([1]) are assumed to obey an expansion in 
terms of the parameter s like 

00 00 
P = i3(0) + J2P(n)£ n , * = *(0) + £ $ (") £ " ' 

71 — 1 n—1 

OO OO 

h A B = «AB + J>j& n , C /A = E C/ (V' ( 2 ) 



n—1 n—1 

where qAB is a unit sphere metric with respect to the 
coordinates x A , and we denote its associated covariant 
derivative by D Al i.e. D A qsc = 0- Note the symmetric 
tensors 7^ contain only two degrees of freedom and arc 
traceless because of the determinant condition on hAB- 
The metric functions have to obey regularity conditions 
at the vertices of the outgoing null cones, because the 
vertices trace the origin of a Fermi normal coordinate 
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system |2l| along the central geodesic. According to the 
vertex lemma in [22j that lists the regularity properties 
of a Bondi Sachs metric at the vertex, we require the 
metric functions to have the following limiting behavior 
at r = 

O(e ):{/3 ,$o} = 0(r 2 ) , (3a) 
0(e n ) : { 7 i"2, £(„),*(„)} = 0(r 2 ) , U A n) = 0(r).(3b) 

Since the metric functions are given by a power series 
in terms of s, any quantity derived from the metric g a b(s) 
is given by a power series in terms of e. For any tensor 
7~(e), we introduce the notation 

oo 

T( £ ) = £f £ ", (4) 

n=0 



where 7" is T(e) evaluated in the background spacetime 
n , 

and 7" is the n perturbation of T{e) with respect to the 
background spacetime. 

With 1Z a b denoting the Ricci tensor, the vacuum Ein- 
stein equations lZ a b = for the line element ([J), can be 
grouped into three sets, a group of three so-called supple- 
mentary equations (1Z UU = and lZ uA = 0), one trivial 
equation (lZ ur = 0) and six main equations consisting of 
four hypersurface equations 

TZ rr {e) = , TZ rA {e) = Q , (5a) 
n ( 2D)(e) := g AB (e)K AB (e) = , (5b) 

and two equations 

Ull(s) := K AB (e) - -g AB (e)1l (2D) (s) = . (6) 

This grouping is given by the Bondi-Sachs lemma [H-Q 
obtained from the twice contracted Bianchi identities: // 
the main equation hold on one null cone and if the op- 
tical expansion of null rausfldi] does not vanishes on the 
cone (i.e. [3 is finite), then the trivial equation is ful- 
filled algebraically and the supplementary equations hold 
provided they hold at on radius r > 0. Therefore, the 
supplementary equations can be seen as constraint equa- 
tions for the metric functions at given radius r > 0. As 
we intent to solve the Einstein equations on the entire 
null cone including its vertex, the regularity conditions 
at the vertices can be used to replace these constraints. 
Although the quasi-spherical approximation being intro- 
duced is completely general, we consider for simplicity 
hereafter only the zeroth and first order terms of in the 
e-expansion of the main equations, whose relevant Ricci 
tensor components are given in App. [XJ 



III. SOLUTION PROCEDURE FOR THE 
BACKGROUND SPACETIME AND LINEAR 
PERTURBATIONS 

A. Background spacetime 

Here we show that the background spacetime must be 
Minkowskian, when a regular vertex and a vacuum space- 
time arc assumed. 

o o 

From lZrr= 0, TZ(2D)~ (App. IA and the regularity 
conditions (|3"a|) it follows that 

Po(x a ) = , $o(x a )=0 . (7) 

Hence the background metric is Minkowskian with re- 
spect to outgoing null coordinates and has the line ele- 
ment 

ds 2 = -du 2 - 2dudr + r 2 q AB {x c )dx A dx B . (8) 

Hereafter, we set /3o = $0 = in all equations and use the 
standard spherical coordinates x A = (9, tfi) to parameter- 
ize the unit sphere metric, i.e. q A B(x A ) = diag(l, sin 2 9). 

B. Linear perturbations 

1. A master equation for vacuum perturbations 

In this section, we derive a covariant differential equa- 
tion for a master function that allows us to determine all 
linear perturbations of the Minkowski background space- 
time. 

From (|A3|) and the regularity conditions (|3bj) it is seen 
that 

[3 {1) (x a )=0, (9) 

which is hereafter imposed in the calculations. Setting 
Tl rA = (from eq. (|A4p) yields 

= (r 4 q AE Uf llr ) r + r 2 q EF D El f Fr . (10) 

From this equation it is clear that if cither 7^ or UA\ 

is known on the interval [0, 00), eq. (|10p can be used to 
solve for the other respective field. Inserting Q and © 

1 

into (|A5|) . and setting TZ[2D)~ gives the equation 

= ^,.-D A D B lA 1) B -^D A (r 4 U A ), r , (11) 

which allows us to determine an intermediate variable $ 
that is related algebraically with by 

<£ = 2r(l + 2* (1) ) . (12) 

The two equations (fTU|) and (fTTj) link the three perturba- 
tion variables. In particular, UA\ and can be calcu- 
lated once 7^ is known or 7^ and $(!) are determined 
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from In what follows we determine and <I>(i 



'(i) 

from 7^3. 

Denning the functional ^{Tab) and £$(DaX e ) act- 
ing on 7^ and arbitrary covectors Xa, respectively, like 



^( 7 <&) := Kr7ilj,r-5(r 2 7il r ) r ,(13a) 
1(D a Xb) := D A XB-\qAB{q EF D E X F ) , (13b) 



allows us to write the evolution equations 1Z A b = ® ( e 1- 
([A6])) briefly as 

= ^(7is)+^^ (qaeDbU^ +2.^ (q AE D B ufy . 

(14) 

Taking the covariant derivative Db of (|10p we can derive 
= r 2 .%(q AE D B U E )rr ) +4r@(q AE D B Uf 1)tr 



+ <%(D B D El V, r ) ■ 

From the following calculation 



(15) 
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d_ 

dr 



a- us 



r-W 



r 2 m , (16) 



we obtain the equation 



(17) 



Defining the auxiliary variable xab '■= r JAB> usm g the 
definition (|13aj) allows us to write (fT7|) as 

,/ 1 







r [XAB.rru " -XAB.rrr 



(DaD 1 



XBE,rr 

(18) 



This a second order differential equation for ipAB 

XAB,rr, »■ e. 







1 



r 4 (2^AB,« - 1pAB,r) 



- 2DaD k iPbe 



q A B[D E D F i> EF 



(19) 



Since a solution of equation (p~9|) is the starting point 
to determine all non-trivial metric fields, we shall call it 
the master equation of linearized vacuum perturbations 
on a null cone. The two tensor field ipAB, which is deter- 
mined by this differential equation, shall be referred to as 
the master function of the linearized vacuum perturba- 
tions. By calculating the linearized Riemann tensor with 
respect to JT]), we find that ipAB determines the compo- 
nents IZrArB via 



n 



rArB — 



-xripAB- 



(21) 



2. Representation of the perturbations and their equations 
in a spin frame work 

In equations. (fTU|) . (fTTj) and (fT^)l the angular deriva- 
tives are covariant derivatives on a unit sphere. In prin- 
ciple, these covariant derivatives could now be expressed 
by the corresponding partial ones utilizing the represen- 
tation of the unit sphere metric in terms of the coor- 
dinates 9 and </>. However, we follow the approach of 
[l5| that became standard when working in numeri- 
cal relativity with the Einstein equations in Bondi-Sachs 
formulation. On the two-surfaces du = = dr (r > 0), 
we introduce a complex dyad q A and its complex con- 
jugated q A to represent the unit sphere metric qAB and 
its corresponding covariant derivative V a- The dyad is 
defined by q A qA — q A u = q A r = and q A q~A = q with 
q > 0. The latter definition covers both commonly used 
normalizations, the traditional one q = 1 [ToL IT5L l26l. \2l\ 
and the numerical one g = 2 @, 0, [28|, which is em- 
ployed in numerical relativity. The unit sphere metric, 
qAB, represented by the dyad is 



qAB = - {qAqs + qAqs 

a V 



(22) 



According to our choice of angular coordinates, x A = 
(9, (f>), q A can be expressed as q A = (g/2) 1 / 2 (l, isin -1 9) 
with i = v 7 — 1. 

Any traceless, symmetric tensor T)rAi...A,) °f rank s on 
the two-surfaces du = = dr (r > 0) can be expressed 
by the dyad q A and a complex scalar field T) like 

V(A 1 ...A S ) = —(v q Al ■ • -1a s +V qA 1 ■ ■ .«a.) , (23) 
where 



Based on our previous definition of the intermediate vari- 
able xab, we find the following equation 



(nAB)-rr = ^ AB 



(20) 



which allows us to determine Jab f rom ipAB- The re- 
maining field perturbations UAs and $(1) are then ob- 
tained by integrating hicrachically (fTO)) and (jTTj) . 



V = V(A 1 ...A s )qA 1 ■■■qA,, 

V = V(A 1 ...A s )q Al ■■■1a, 



(24a) 
(24b) 



When the spacelike vectors $ie(q A ) and ^sm(q A ) are ro- 
tated in their complex plane by a real angle the quan- 
tity 7] transforms under this rotation as 7/ = e^rj and 
it is said to have the spin weight s [HI [HI • For the co- 
variant derivatives Dbtj^Ai...a,) with respect to the unit 
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sphere, we define the eth and eth-bar operator as 
[2 



9?7 := 
5r/ := 

which correspond the following derivatives of 77 in terms 
of the coordinates x A 



q Al ...q^q a D B j 1{Al ... As) (25a) 
q Al ...q A *q B D B r] (Al ... As) , (25b) 



r) = (sin s 9) 
1 



sin s 6> 



8 i 8 
89 sin 9 8(f) 
8 id 



V 



(26a) 



89 sin ( 



sin s 6» / 
(r]sm s e) . (26b) 



Based on (|25|) . we find the commutator of the eth and 
eth-bar operator as 



[3, S]t7 = 2sK.r] , 



(27) 



where K, — 1 is the Gaussian curvature of the two- 
surfaces du = = dr, (r > 0) with respect to the 
unit sphere metric, which is calculated from K, := 
(l/q 2 )q A q B q E q F n ABEF (q CD ). 

To write (fTU|) . (fTTj) and (jT^J) in a spin representation, 
we define the spin weighted quantities 



ip := q A q B i' AB , 
W := g^E/ft , 



(28a) 
(28b) 
(28c) 



where U has spin- weight 1 and ^ and J have spin- weight 
2 [4l[ . According to ([24aJ , the metric perturbations U A 



and ~f AB can be found from U and as 
1 



(1) 

lAB - g2 



Uq A +Uq A 



JqAiB + JqAqB 



(1) 

(29a) 
(29b) 



Multiplying (fT9|) with yields the spin representation 
of the master equation 







1 



4U 



(30a) 



From the definition of iJj ab wc deduce an equation that 
determines J from ijj 



(rj) =1/J . 

\ / .rr 



(30b) 



Multiplying (|T0|) with q A and using (|22[) allows us to 
write (fTOl) as 



= r 4 W, r + 



(30c) 



The spin representation of cq. can be found as 




1 



r 2 V2q 



m 



(30d) 



where 4> has a spin weight zero, because it is a tensor of 
rank zero. Although it is not required hereafter, we give 
for completeness the evolution equation (|14p in terms of 
the spin weighted variables 







2(rJ), 



r 2 J, r 



\2q) 1 / 2 dU 



(31) 



In equations (|30cp , ()30d[) . and (|3"Tj) it is seen that the 9 
and 9 operators raise and lower, respectively, the spin 
weight when applied to a spin weighted quantity, since 
the overall spin weight in these equations must coincide 
with the spin weight of the variable that does not carry 
an 3 or 9 operator. 

Equations ([3U]) are the equations that need to be solved 
for the perturbations in the spin representation. Equa- 
tion (|30a|) is a wave equation for the spin-2 field ip and 
we show in App. [Blhow it is related with the flat -space 
wave equation of a spin-0 field. The field ip is, in fact, 
the linearized Weyl scalar 2r\E'o [l?}. Since a solution of 
(|30a| deter mines all linear perturbations with respect to 
a Minkowskian background, the perturbations describe 
the propagation of spin-2 waves on a Minkowskian back- 
ground spacetime. 



Decomposition of the perturbations into spin-weighted 
harmonics 



A standard approach to solve the linearized Einstein 
equations is to decouple the angular dependence from 
the equations by expressing the perturbations in terms 
of a spectral basis depending on the angular coordinates, 
only, and where the coefficients of such spectral series 
depend on all other coordinates but the angles. By in- 
serting such spectral series' into the equations, one then 
is able to derive differential equations for the coefficients 
of these series. Hereafter we follow this approach. 

Since in eqs. (|30a|) . (|30c|) and (|30dj) the angular deriva- 
tives are expressed in terms of the 9-opcrator, we de- 
compose "0, J, U and $ with respect to a basis of spin 
weighted harmonics s Yi m which are eigenfunctions of the 
operator 99. Let Yi m (x A ) be the conventional spher- 
ical harmonics [29], then the spin-s weighted spherical 
harmonics, s Yi m (x A ), are derived from the spherical har- 
monics, Yi m (x A ), like [2f| 



s^lra 
sYltn 



k(l, s)<5 s Y lm for 
Y [m for s = 

(-)l s lfc(/, \s\)Q lsl Y h 



s > , (32a) 
, (32b) 

, for s < , (32c) 
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where k(l,s) := [(I - s)\/(l + s)!] 1 / 2 and s Y lm = for 
\s\ > I. The follo wing: properties of the 3-operator and 
the s Yi m are used [421 ]: 



sY Ijn 



( _ r+s _ s Y K _ m) , (33a) 

+ y/(l-s)(l + s + l) a+1 Yi m , (33b) 

-V(i + s)(Z-s + l) s _i^ m , (33c) 

-(Z-s)(Z + s + l) s Yi m , (33d) 



where (|33d[) shows that s Yz m are the eigcnfunctions of 
the operator 33 [43[. In particular, if s = the operator 
33 corresponds to the angular momentum operator since 
33Y; TO = —1(1 + l)Yi m . We make the following ansatz for 
ip, J, and U 



iP(x a ) 
U{x a ) 



Y Y ^i™{u,r) 2 Y lm (x A ) , (34a) 

(=2 m=-l 
oo / 

Y Y Jim{u,r) 2 Y lm {x A ) , (34b) 

1=2 m=—t 
oo I 

Y Y Ui m {u,r)iY lm (x A ) . (34c) 



l—l m——l 



The perturbation variable $ is a spin-0 field, thus it 
would be the most natural to express it in a oYi m ba- 
sis. However, an inspection of (|30dj) while using ([33]) 

shows that the terms containing 3 J + B 2 J and WA + WA 
have an angular behavior like 



oYlm + (— } 0^(-m) 



By defining the following spin-0 harmonic 



Zlm(x A ) 



oY lm {x A ) + {-) m Y l{ _ m) {x A ) (35) 



and making for 4> the ansatz 



$(x a ):=Y Y ®lm(u,r)Z lm {x A ) , (36) 



l— 1 m— — l 



allows us to decouple the angular derivatives in (|30d|) 
with the angular part of $. 

Inserting (JMJ) and (J35J> into ([HI), and using (|33p yields 
an hierachical set of differential equations coupling the 



coefficients ipi m , Jj m , C// m , and 
1 



rJi 



= — |r 4 (2V>i TO ,„ - 1pl m ,r) 

+(i + 2)(l-l)^ m , 

V'im , 



, /(Z + 2KZ-1) 
V 2 <? 

l)Z(Z + l)(Z + 2) 
9 



(37a) 
(37b) 

(37c) 



(37d) 



IV. BOUNDARY CONDITIONS AT THE 
VERTEX 

The metric perturbations variables, J ', IA, and $, are 
functions on null cones u = const with their vertices 
traced by the central geodesic of a spherically symmetric 
background spacetime. The null cones, however, are not 
diffcrcntiable at their vertices, consequently the metric 
perturbation variables are also not differentiable there. 
As the perturbation J defines the complex master func- 
tion ip, the variable ip is a priori also not differentiable 
at the vertex. In this section we follow the procedure 
as in [22j to discuss the boundary conditions for ip and 
the metric variables at the vertex. Thereby, we assume 
a Minkowskian observer and the existence of an convex 
normal neighborhood along the central geodesic of the 
spherically symmetric background spacetime. 

Since the master function ip is determined from per- 
turbation 7^ like ip — (l/2)(rj^q A q B ) rr and since 

7ab = 0(r 2 ) at r = 0, we conclude the limiting radial 
behavior of ip to be 



ip{x a ) = 0{r) 



(38) 



To find the exact behavior of ip near r = 0, we assume 
that ip obeys, at r = 0, an infinite power series expansion 
for ip in terms of r like 



iP(x a ) 



n=l 



ip n (u,x A )r n 



(39) 



where the coefficient functions ip n (u, x A ) are assumed 
smooth functions for all values of u and x A . Accord- 
ing to the vertex lemma of [22|], the radial coefficients of 
(|39| must obey the field equations. Inserting (|39|) into 
(|50a| yields 

OO s 

= -r33V-i-^|(n + l)(n + 4)^„ +1 

n — 1 ^ 



+ (33-4)V„+i-2(n + 4)^„ jU U' 



(40) 
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For this equation to be fulfilled, we deduce the conditions 

o = gg^i , (4ia) 

= (n+l)(n + 4)^ n+1 + (gg-4)^ n+1 

-2(n + 4)^ n ,„ , (n>l) . (41b) 

To decouple the angular and time dependence, we make 
the ansatz WJI 



-l i 



l(j n (u,X A ) = ^ E C l-m.n{u) 2 Yi m {x A ) 



(42) 



1=2 m=-l 

Inserting (|42j) into (|41a|) yields 

2 



o = ggvi= E (-5)(2-2)c 2 . m .i(«) 2 y (m (^) , (43) 



from which we conclude that the functions c n +2.m.n+i(u), 
(n > 1) are 2n + 5 arbitrary functions for the to— modes 
of 2^(n+2)m and the coefficient functions (u) are 

for (71 > 1), (2 < I < n+l), \m\ < I given by the recursive 
relation 



Cl.m.n+l(u) = a Ln —Cl. m . n (u) 

du 



where 



a-i.n '■- 



2(n + 4) 



(46) 



(47) 



(l + n + 3){l -n-2) 
Defining the arbitrary functions c n +2.m.n+i(u) as 

C n +2.m.n+l{u) ~ C^ 2 n+1 {u) , (jl > 1) , (48) 



showing that c 2 . m .i{u) are five functions for the 
to— modes of the I = 2 spin-2 harmonic, 2^2m, that can 
be chosen arbitrarily, and which we write as 



C2.m.l{u) := C^™\u) 



(44) 



Inserting (|4"2"]) into (|41bj) gives, after little algebra, for 
n > 1 



n+2 

= c n+2.m.n+l 
m=-(n+2) 
n+l I 



4-4 



2^(11+2)171 



, Cl.m.n 

du 



+E E - 2 (™+ 4 ) 

1=2 m=-l 

+ci. m .n+i [(n + l)(n + 4) - (/ - 2)(Z + 3) - 4] J 



(45) 
I 



and writing out the recursive series (|46[) for the first val- 
ues of I and n allows us to deduce for (n > 1), (2 < I < n), 
I ml < I an explicit form of the recursive series (|46p 



jn-l+l 

ci. m . n {u) = ( I] a l . k ) 1 ^ m C\jl 1 {u) . (49) 



k=l-l 



The product term involving the ai.k in (|49p can be sim- 
plified further using (|4"7|) and properties of the factorial 
and Gamma function 



ra-1 

ai.k 

k=l-l 



(2Z + 2)!(n + 3)! 



(1 + 1) (n + Z + 2)!(n-Z + l)! 



(50) 



Using (32), gl]), (gll), g6]), (J48]), dMl) and (HOj allows us 
to find the explicit dependence of ip from the arbitrary 
functions C^ +1 , that is 



00 n+l 



V>0O = E E [cl%M2Y {n+1)m (x A )]r n 

(2Z + 2)!(n + 3)! 



n=l m=-(n+l) 
oo n I 

EE E 

n=2 £=2 m=-l 



(Z + l)(Z + 2)!(n + Z + 2)!(n-Z + l)! 



(51) 



Equation (|51[) shows the boundary conditions of ip in 
terms of a power series in r at the vertex and functional 
dependence of ip from data - the functions C^ +1 (u) - 
that are given as free functions along the central geodesic. 
These free functions are the spin-2 multipoles of the com- 



plex master function, i.e. the Weyl scalar 2r\E'o, which 
are calculated with respect to a Fermi observer follow- 
ing the central geodesic. Given ijj we find the boundary 
conditions for the perturbation J, U, and <J>(]\ as 



8 



OO I 



j^ a ) = E E [##(«) 2 w**) 



oo n— 1 / 



1=2 m=-l 



n=3 1=2 m=—l 



E b ^ o 



jn—l 



— J&°(u) 2^00 



(7l/ 



(52a) 



W(:r Q ) 



2. 2..[ ( i + 2 ) 72I 7 " 



£=2 m=-Z 

oo n— 1 i 



j-i 



EEEtt 



i6(«i0 \/( z - 2 )(' + 3 ) 



r rin-l 



n=3 1=2 m= — l 



(l-n)(n + 2) 



2g 



(52b) 



OO I 



= E E [ 



vg-iw +w+ 2) 4:; i) («) ^™(^ A ) 



Z=2 m=-l 

OO 71— 1 i 

+EE E 

ri=3 i=2 m= — l 



Aq(l + 1) 

^(l-l)l(l + l)(l + 2)b(n, I) 
iq(n + 1) 



(52c) 



where we defined 
6(n,0 



n(n + l) ' 

„„_!_, Z(n + 2)(2; + 2)!(n 



D! 



(l + 2)\(n + l + l)\{n-l)\ 



(53) 
(54) 



It is seen in (|52ap that for static spacctimes the free data 
Ii™\u) of the I— multipole the along the central geodesic 
correspond with power r l , and for dynamical spacctimes 
the n th time derivative of the I— multipole of the free data 
are found at the power r n+l . Moreover, formula (|52a|) 
agrees in axial symmetry (i.e. m =0) with expressions 
given in [22]]. 



V. SOLUTIONS APPROACHES TO THE 
MASTER EQUATION 

In this section, we present two different approaches to 
solve the differential equation for the spectral coefficients 
of the complex master function ip. In the first approach 
(Sec. IV A[) . we make a 'standard' ansatz to solve the dif- 
ferential equation, whereas in the second approach (Sec. 
IV Bl) wc impose an ansatz based on the characteristic 
nature of the partial differential equation. In both ap- 
proaches we look for solutions which respect the regular- 
ity conditions of the previous sections and we also provide 
the initial data for the free functions C^?\ n (u) . 

A. First Approach 



To solve the differential equation (|37a[) for the coeffi- 
cients of ipim, we make a 'standard' ansatz of separation 



of variables 

^i m (u,r) = e T ^R lm (r) 
Inserting ([55]) into (|37a|) yields 



dT 
du 



4r 



dRi 7r 
dr 



1(1 + 1) -2 



R 



8rR lm + 2r 2 



dRi„ 
dr 



(55) 



(56) 



Since the right hand side of (|56j) is independent of u, it 
can be treated as a constant in the u— integration that we 
set to be dT / du = B = const. This implies the solution 
T(u) = Bu + uq, where the constant uq is determined 
by the initial conditions, and without loss of generality 
we may set uq = 0. Inserting dT/du = B into (|56|) and 
rearranging gives 



= r 



d 2 Ri m 
dr 2 



r 2 B)^^-(8rB+l(l+l)-2)Ri m 



dr 



(57) 



To find a solution for (|57j) . we make the ansatz Ri m (r) 
r k e ar Ai m (r) and obtain 







^ d 2 Ai m (r) 

dr 2 
a(a - 2B)r 



2 (a - B)r 2 + 2(k + 2)r 
(4a + 2ka - 2Bk - 



dA hn (r) 
dr 

A im (r) 
(58) 



-l(l + l) + 2-3k-k 2 A lm (r) 



If we choose k = —1, a = B, and make in (|58| the 
variable transformation z = Br, (|58|) can be cast into a 
inhomogencous Bcssel type differential equation 



> d 2 A h 
dz 2 



dAv 



dz 



z 2 + l(l + l) 



A, 



AzAi 



(59) 
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whore the inhomogenity is given by the right hand side of 
([59| . The homogeneous counterpart to |59|) is the Bcsscl 
differential equation 

,d 2 A lm (z) dA lm (z) 



~d? 



+ z- 



dz 



z A + 1(1 + 1) A lm (z) = Q 



(60) 

whose solutions arc the modified spherical Bessel func- 
tions, ii(z) and ki(z), of first and second kind, respec- 
tively, 



ii(z) 



h+dz) , h(z) 



K l+h (z) , (61) 



where I l+ i (z) and K l+ i (z) are the modified Bessel func- 
tions of first and second kind [3(| P- 437]. The regularity 
conditions (f5Tj) for ip require that Ri m (r) = 0(r' _1 ) for 
r w and I > 2, consequently Ai m (z) must behave as 
0(z ) for z « 0. Since the functions i;(z) and fc;(z) be- 
have as O(0 Z ) and 0(z - ') for z sa 0, respectively, the 
modified spherical Bcsscl functions of the second kind 
must be ruled in the construction of a regular solution to 
(l59l). To solve (l59l). we make the ansatz 



Al m {z) 



2 

£ 

fc=0 



where are the coefficients , and , are determined 
by the boundary conditions at z = and by inserting (|62|) 
into . Inserting (|6"2"j) into (JSHJ) allows us to deduce the 
solution 



Alm{z) = \Z 2 + Z)k{z) 

+./f™ ) \z 2 + (l-l)z + 



1(1 -I) 



k-x(z) 



(63) 

Thus a master function ip with the following spectral co- 
efficients 



4>im(u,r) 



f: 



(Br) 2 + 
1(1-1) 



Br 



ii-i(Br) 



ii(Br)\ (64) 



(Br) 2 - (I - l)(Br) + 



respects the regularity condition (|51|) at the vertex. In 
particular, if we choose in ([STjl the free initial data to be 



c^Au) = ft l) 



1 



1(1-1) 



,fltl Tjl 



.(2Z + 1)!! 2(2/ + 3)!! 

where /!! is the double factorial, then the free data ([63)1 
describe the solution (|64)) in a neighborhood of r = 0. 

It can be easily seen in (|64|) that the spectral coeffi- 
cients ipi m diverge exponentially as r tends to infinity. 
Consequently the spectral coefficients of linear pertur- 
bations that are calculated from (|64|) will also diverge 
exponentially and hence the solution which we found for 
the master functions will not give to an asymptotically 
flat spacetime. 



B. Second Approach 

In this section, we present an approach yielding a solu- 
tion of the master equation that is regular at the vertex 
and finite, when r tends to infinity. 

Defining -0; m = (1/r 4 ) J r 4 Pi m dr shows that (|37a|) is 
an intcgro-diffcrential transport equation 







Pi, 



1 



1(1 + 1) -6 
2^6 



r 4 P lm dr .(66) 



It can be seen that the integral part of (|66[) vanishes for 
/ = 2 implying the surfaces u + 2r = const to be the (in- 
going) characteristic surfaces of this integro-differcntial 
transport equation. 

Since eq. (|37a[) is another representation of the trans- 
port equation (|66|) and based on the fact that u + 2r is 
a characteristic surface for the lowest (1 = 2) multipole 
of ip, we make the following ansatz to solve (|37a|) 



iplm = Ci m u l ri(u + 2r) k 



(67) 



where the exponents i, j, k need to be determined by 
inserting (|67| in (|37ap and Ci, n £ R are arbitrary con- 
stants. The calculation yields the following possibilities 
for the exponents for (I < 2) 



(1 + 2), -(1 + 2), (1-2) 
-(l-l), (l-l), (1 + 3) 



[«, j, k] 
[i, J: k] 



By inserting these values into (|67|) and expanding the 
thus obtained function at r = 0, it is seen that the first 
possibility, (|68a|) . for the exponents gives rise to a sin- 
gular behavior of the coefficients V'Zm in (|67|) at r = 0, 
whereas the second one complies with the regularity re- 
quirement (|5ip . Hence, ip with following spectral coeffi- 
cients 



1plm(u,r) = C h 



2r) 



1+3 



(69) 



yields for u > a complex master function ip that is 
regular at the vertex at r = and also finite when r — > 00 . 
The initial data along the central geodesic for u > are 



°Z+l.i 



Clr, 



u 2(l+2) 



I > 1 



< I 



(70) 



B , (65) VI. 



ASYMPTOTICALLY FLAT AND REGULAR 
SOLUTION FOR THE VACUUM 
PERTURBATIONS 



Given the spectral coefficients (|69|) for the complex 
master function and the boundary conditions for the 
perturbations (|3b)l . we now integrate (|37b|) , (|37c[> . and 
(|37d|) for the spectral coefficients of the perturbations. 
Together with f3(i) = and (fT2^) . our solution of the 
linearized perturbations J , U, and $0 i n term of spin- 
weighted harmonics are 
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oo I 



U{x a ) 

<£ ( i)(z a ) 



\- Chn r l [(l + 2)u + 4r] A 

^ ^ Z(l +!)(/ + 2) u l + 2 (u + 2r) l +i 2 lm[ ' 



1=2 m=-l 

oo I 



J_f^ C lm r l ^{lu + 2r) A 

V 7 ^ h m t±z Kl + 1) V^i - 1^ + 2) " i+3 (« + 2r) ; 1 ' 



oo i 

i-iy y 

9 An ^ ^ 



C h 



l - j 2n f- l V(i-mi + m + 
i — 



) tr+^u + 2r)' 1 



(71a) 
(71b) 
(71c) 



which is well defined for u > 0, and where the coefficients 
C/ m are arbitrary constants. The solution above de- 
scribes for u > spin-2 fields propagating on a Minkowski 
background spacetime, and we shall therefore refer to it 
hereafter as SPIN-2. 



A. Asymptotic properties 

1. Conformed compactification, frame at null infinity 

The study of the asymptotic behavior of solutions of 
the Einstein field equations |45| is most elegantly done 
by using Penrose's definition of asymptotically flat space- 
time [23l ] . The idea of this compactification is that, aside 
from a "physical manifold" ^ with a "physical metric" 
g a b, there exists a positive function £, which decreases 
along all complete null geodesies, approaching to zero as 
their affine parameter goes to infinity. Thereby a "non- 
physical metric" g a b = £ 2 f/ab can be extended smoothly 
to a larger, compactificd manifold M = ^ U 31 1. 
The boundary J? := dAi is called null infinity and one 
has I = and V a £ ^ on Points on J? in the man- 
ifold ^# correspond to "points at infinity" for radiative 



fields in the physical manifold. It can be shown that the 
boundary <f is a null hypcrsurface, and that the Weyl 
tensor, C a bcd- behaves as 0(£) in the neighborhood of Jfi 

To find a compactified metric g a b for a given Bondi- 
Sachs metric of a physical spacetime, we have to find 
a conformal factor I such that it has the properties at 
as stated above. Assuming g a b is smooth in we 
first define coordinates x a as a function of the Bondi 
Sachs coordinates x a such that 'points at r — oo' in the 
Bondi-Sachs coordinates are located at finite values in x a 
coordinates. Second we calculate an intermediate metric, 
g* b say, by transforming the Bondi-Sachs metric to the 
conformal coordinates x a and finally we find a conformal 
factor i making £ 2 g* b finite at the location of null infin- 
ity. The coordinates i a are called conformal Bondi-Sachs 
coordinates. 

The coordinate x° := u and the coordinates x A := 
x A are defined as their counterpart in the Bondi-Sachs 
coordinates. The coordinate x 1 := x ranges over the 



interval [0, a], a > 0, and is connected to the physical 
radial coordinate r via a strictly monotonous positive 
function r = r(x) having the further properties r(0) = 
and ]hn^ a r(x) = oo. The points at null infinity J? are 
located at x a = (ii, a, x A ) where u and x A are arbitrary. 
The requirement Ymix->a r{x) = oo and the monotonicity 
imply that r(x) has a singularity at x = a which is a pole. 
Hence the most general form of r{x) may be assumed to 
be of the form 



r(x) 



R(x) 



(a - x) r 



(72) 



where m £ {1, 2, ...} is the order of the pole of r(x) 
and R(x) is a strictly monotonous positive function with 
i?(0) = and R(a) ^ 0. If one now transforms the 
metric components of the line element ((T|) to confor- 
mal coordinates x a and assumes that metric functions 
fiAB, /3, U A , and <& have only poles at x = a, it can be 
seen that the thus obtained metric g* b has a poles of the 
order 2m and m + 1. Assuming a conformal factor like 



I = (a — x) k / 2 , k = {1, 2, ...} and requiring g a b 
to be finite as x — > a implies 



k 



l. 



9ab 



(73) 



Thus the function t[x) = a — x vanishes for x = a and 
its derivative does not vanish at x = a. Therefore the 



surface with the coordinate values 



l ), with 



arbitrary values of u, and x , corresponds to null infinity. 

The conformal metric g a b has the non-trivial covariant 
components 

gau = -(a- x) 2 e 2<s>+4 P + R 2 (x)h A6 U A U § (74a) 



, „ . dR 
(a-x)—+R(x) 
ax 



9uA = -R 2 (x)h AB U B 
9ab = R 2 (x)}iab ■ 

and the nonzero contravariant components 

-l 



fJ 



, „ , (ii? „ , „ . 
(a-x)— +R(x) 
ax 



-2ri 



(74b) 

(74c) 
(74d) 



(75a) 
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(a — x) re 



2„2* 



(a-x)§ + R(x) ' 
l> A [(a-x)^| + i?(x) 



-2/3 



5 AS = 



1 



U 2 (x) 



.IB 



(75b) 



(75c) 
(75d) 



In the quasi-spherical approximation of SPIN-2, the 
conformal metric g ab is expanded in terms of the 
s— parameter as 



o 1 
9ab= 9ab + £ 9ab i 



where g ab has the non-vanishing components 
= -(a-x) 2 , 



" 

9uu 


9ux 




(a — x)—— + R(x\ 
ax 



9ab = R (x)qAB 



The non-zero components of g ab are 



9uu = -2(a-x) 2 $ (1 



l 

9ux 

1 

9uA = 



(1) ' 

, „ . dR , „ , 
(a-x)— + R(x) 
ax 

^R 2 (x)q AB U^ ) , 



9AB = &@YfAB > 



(76) 

(77a) 
(77b) 

(77c) 

(78a) 
(78b) 

(78c) 
(78d) 



where U A ^ and j AB are calculated from hi and J using 



If we calculate g ab at null infinity, i.e. set x = a, and 
define new coordinates x a as 

x := u = u , x := x = — K(a)x , x := x = x , 

(79) 

then a coordinate transformation of ([77j) at from x a 
to x a brings (|77|) to a conformal Minkowski metric g ab 
with the non-zero components 



fjo 



(x Q ) = 1 , 9ab (x a ) = OAS , 



(80) 



which corresponds to the metric of an inertial observer 
in the (compactified) conformal spacetimc (Ai, gab) [II- 
The inertial conformal frame is the frame at null infinity 
where one can uniquely define the asymptotic properties 
of an asymptotically flat spacetime, because most gen- 
erally a frame at null infinity is not a Minkowskian one 
[32| . The Bondi-Sachs variables have in an inertial frame 
at null infinity the values 

h AB \s = q AB , = U A \s=P\s = &\s , (81) 
which means that the perturbations /3(i), U, J , and tfrm 
have the values 



= t 7U=Z?U = ) S ( i ) U = $ (1) 



(82) 



Inserting (|72|) into ([71) and expanding the thus obtained 
expression at x = a yields 



J(i a ) 
U(x a ) 

*(i)(i a ) 



E E o/-x 



Cl m2 Ylm{x A ) 



1=2 m=-l 



2l-W+H(l + l)(l + 2) 



O 



(a-x) 



E E /-i 



C im iY lm (x A ) 



=2 m=-i 
00 £ 

EE 

1=2 m=- 



' 2 J - 1 u'+ 3 1(1 + 1)^(1-1X1 + 2) 



i?(a)C im Z im (x A ) 



, (a - x)2 J + 1 u'+VC- +W+ 2) 



(a-x) 2 
+ o[(a-s)° 



(83a) 
(83b) 
(83c) 



Expressions ([83") show that J ,U 1 and $m evaluated at 
exhibit non-trivial values which indicate that SPIN- 
2 does not result in an inertial conformal frame at null 
infinity after applying Penrose compactification. Note 
although (|83c[) diverges at J , the corresponding metric 
components gnu and g xx are finite at ^ . 



2. On the determination of ^4 in a conformal frame 



To extract unambiguously physical information from 
SPIN-2 at null infinity, we have to find, in principle, a 
coordinate transformation that casts the metrics ([77)) and 
(|78p into a conformal Bondi frame at null infinity. In this 
section, we determine the We yl s calar "J 4 [13, in a 
Bondi frame at J? by following [24j . In this procedure, we 
do not look for the exact coordinate transformation from 
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the conformal frame to a Bondi frame, but start with an 
expression of ^4 in a Bondi frame and then transform it 
into an arbitrary conformal frame at J? . We be gin with 
a summary and motivation of the basic ideas of [24| and 
give the expression of ^4 in a Bondi frame at J* for the 
solution SPIN-2 in the end. 

From the conformal metric g a b = £ 2 gab in and the 
vanishing Ricci tensor, 7Z a b(g a b), in the physical manifold 
the following equation can be derived 



= l 2 K ab + 2e\7 a \7 b e+gat 



V C V C ^-3(V C £)(V C £) 



(84) 



Taking the trace of (|84|) and subsequently the limit (£ — > 
0) of this trace, shows that the surface with the normal 
vector V c £ is a null hypersurface and that (V c i?)(V c £) — 
0(1) . From (|84|) we can derive two equations 



= en + 6V a v a e- 
Hn ab -jg ab n 
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(85a) 



V a V a l--g ab (X7 c VJ) 



(85b) 



where 1Z is the Ricciscalar with respect g ab and we used 
(|85a| to eliminate the term (V C £)(V C ^) in flM} to obtain 
(|85b[) . Evaluating ([55]) at J?, allows us to define two 
fields that arc intrinsic to J* 



where 



6 := V a V a ^ 



S ab := V a V Q £- -? ab (V c V c i 



(86a) 
(86b) 



(87) 



and which have been first found in Q , and are crucial in 
the calculation of ^4 in (24|. Using i = (a — x), we find 
as 







1 



-99 



and from (fTH jl . lj8Tjt . and V a V b £ = f £ ab follows 



* x a. - 2. A J_ /f /_\ 

uu — 9 9ua,u 2^ 9uu,a ^9uu™ 3 

.9 5u[^,a] 2^ 9ux,u 



1 1 / ^xu p, 1 %xB * 

— .9 .9u[b : A] + 2V 9 9uA,it-\-g 9AB.ii 



~xb ' 

L 

1 



|&L40 



(89a) 
(89b) 

I 

(89c) 
(89d) 



= 3 9u(x,A) + n5 3AB,i , 



(89e) 



Z-'AB = 9 9u 



(A,B) + 9 9C(A,B) - 2 



1 



|3AB0 



(89f) 



The calculation of ^4 in [24j involves three different met- 
rics: 

1. the Bondi-Sachs metric g ab of the physical space- 
time, 

2. the conformal metric 7} ab that maps null infinity to 
a finite value in the conformal space time, 

3. and the conformal inertial metric g ab . 

As the conformal metric g ab is related to the Bondi-Sachs 
metric with the conformal factor t, the conformal Bondi 
metric g ab is related to the Bondi-Sachs metric with an- 
other conformal factor like g ab = Q 2 g ab - As in [24l.[33|. 
we set 



Cl(u, x, x A ) := £(x)lo(u, x a ) 



(90) 



This choice for the conformal factor VL (together with the 
particular choice of the function r(x) and £(x)) has the 
advantage that the compactification in radial direction is 
decoupled from the angular behavior of the metric at ^ . 
The definition (|90[) relates the conformal metric g ab and 
the conformal inertial metric g ab like g ab = w 2 g a b- 

Suppose x a are inertial coordinates as in the previ- 
ous section such that the metric g ab at <f has the form 
as in (|80| . In this coordinate system we choose two 
real null vectors n a := g ah V b {l\^ and l a := g ab \7 b u\j: 
and a complex null vector Q a that obey l a n a = — 1, 

Q a Q a = q , whereas all other scalar product vanish 
between them. These null vectors define a null tetrad 
z?b) : ~ (^°> ™ a > Q a > *3 a ) allowing us to write the inertial 
conformal metric at J" as 



(Jab = —iafhb - nj b + - [QaQ b + Q a Qb 



(91) 



With this null tetrad, the Weyl scalar ^4 is given by a 
contraction of the Weyl tensor [24], HH 



*4 



— lim 

q n^o 



l a Q b n c Q d C a 



bed 



O 

LP) 



(92) 



which corresponds to — (l/g)^ 4 in the standard 

Newman-Penrose notation [l7|, where ^4 is the Weyl 
scalar as defined in [23| . 

The above defined null tetrad is not completely fixed 
[35L |36| as one still has still the following three freedoms in 
the tetrad representation of the metric (|9ip : (i) Lorentz 
transformations with the boost factor a and spatial ro- 
tations around the angle •& (where a and 1? are real func- 
tions), i.e. 



1 



a 



e M Q a 



(93a) 
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(ii) null rotations around n a with a complex function k 



l a = l a + nQ a + kQ 
Q a ' = Q a + \h a , 



n a , (94a) 
(94b) 



and (iii) null rotations around l a with a complex function 
A 



l a = l a , n a = ? 

T' = Q a + AT" . 



\6 a + \6 a 



(95a) 
(95b) 



where the relation C a bcdh a M b n v h d = was used. The 
main result of [24| and the fourth ingredient for the de- 
termination of ^4 is that (|101[) can be expressed by the 
vector field £ a & like 



1 1 

q a) 3 



n a M B M c (V a E BC -V B E aC 



(102) 



This equation has an advantage to (|101[k it is easier to 
determine from the metric at than the (rather tedious) 
calculation of the contractions of the Weyl tensor. 



Since the Weyl scalar (f9"2")) is invariant under null rotation 
around h a on J" ', null rotations around n a can be used 
to calculate "J 4 in another frame at J? . 

Babiuc et al. [24| use four ingredients to find the rela- 
tion between ^4 in an inertial conformal frame and an 
arbitrary conformal frame on <? . The first ingredient is 
that the inertial conformal metric g a b and the conformal 
metric g a b are related via g a b = u 2 9ab- This implies that 
the Weyl tensor transforms between the both frames as 

12 



Cabcd — & C a 



bed ■ 



(96) 



At ^ , the conformal 2-metric gAB is subject to the con- 
straint gAB = QJfi jQABi which yields the following el- 
liptic equation [33j 



3. Calculation of * 4 for SPIN-2 

We find the Weyl scalar * 4 for SPIN-2 by deriving 
first the corresponding expression to (| 102[) in a quasi- 
spherical expansion while assuming the following limiting 

non-trivial values of g^b anc ^ 9ab 



ffia U = -R( a ) , 9ab\s = R\a)qAB , 



I „ - 



(103a) 



5fiiU = 

y AB 



-R(a) , g^W 



-J? 2 (a)g AB C/ ( ^m03b) 
(103c) 



To determine the conformal factor uj at ^ , we consider 
the quasi-spherical expansion 



(g AB ) =2[u* + g AB V A y B log w 



(97) 



1 

w(e)= w + CJ £ 



(104) 



allowing one to calculate the conformal factor to from 
the curvature scalar 1Z of the surfaces u = const on . 
The second ingredient concerns the expression of the null 
vector n a in the conformal frame at , i.e. 



n a U = -n a 

UJ 



(98) 



where n a is given by 



= g ab v b i , = -g ax 



g ux \ 



0. 



-9 



(99) 

Equation (|98| shows that the transformation between 
the inertial conformal frame and an arbitrary conformal 
frame corresponds to a Lorentz transformation (boost) 
from one frame to the other. The third ingredient is 
to use (I95b[l to transform Q a to an arbitrary conformal 
frame 14611 . i.e. 



Q" 



u> 



(100) 



where M a \y := (0, 0, F A /R(a)) such that M a M a = q 
and M a M a = at J. Combining (92]), ®, dSHJ), and 
(|100p , gives the Weyl scalar "J 4 in an arbitrary conformal 
frame at J? 



1 1 

q a) 3 



lim 



n a M b h d M d C abc d 



(101) 



Calculating zero order terms in e of <?ab(£) = qAB/& 2 (s) 

1 
at J* yields Cj= 1/R(a). To determine ui, we calculate 



the 0(e)— contribution of (|97|) while using d)= 1/R(a) 
which gives us the equation 



2R(a 

or at J? in terms of the 9— operator 



1 



D A D A oj +2 uj , (105) 



AqR(a) 



_ 1 1 

99 cj +2 w 



(106) 



Since Cj\ is a real scalar field it has spin weight zero, 
therefore we assume for uj\ an expansion in terms of Z; m 
like 



00 m 



(107) 



=-1 



Inserting (|T07l) and (|83a|) into (fT06|) while using (|33p and 
([35jl implies the spectral coefficients w Zm 



o/ m (ii) 



; K-D 



£ 10 = c2, u 







(108) 



2l+3fiJ+2 g j 2 _ ^ + i)]^(z + + 2) 
(Z > 1) and \m\ < I . (109) 
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The quasi-spherical approximation of null vector h a can 
be found using (|99j) . i.e. 



n jr = 



-(>- 

(a) V u 



i?(a) 



(110) 



To find M at ^ , we use its quasi-spherical expansion 
.1/0 



M 1 



e F 



(111) 



its normalization M a M a - q = M a M a = 0, the quasi- 
spherical expansion for the metric (|103[) . and the rela- 
tions (|22[) . and (|29b[) which allow us to deduce 



F A = q A 



F A - -f I* 
2g 



(112) 



An inspection of ([89]) at ^ shows that E Q b is O(e) at 
J ' . Therefore, ^4 is of O(e) at J' and only the O(e ) 
parts of w, n a , and M A must be taken into account in 
its calculation. Since we have n x | > =0 and because E a b 
is of 0(e) at ,f, only the following covariant derivatives 
of E a & must be considered at 



^C^uB 



= 



R{a)e (113) 
(114) 



With these covariant derivatives we find ^4 in the quasi- 
spherical expansion as 



— <5U : u + — S,uu 



R(a)e , (115) 



which corrects equation (3.54) in [2J|. Using (|83ajl . ()83bl) . 
and ()115p for SPIN-2 yields a simple expression for \I/4 
at null infinity 



1=2 m=-l 



(I + 3)Cl m 

2 l u'+ 4 l(l + l) 



iY lm {x A ) (116) 



VII. SUMMARY 

We discussed a linearized vacuum solution of the Einstein 
equations in the Bondi-Sachs formulation of General Rel- 
ativity. Assuming that the metric obeys regularity con- 
ditions along the central geodesic tracing the vertices of 
the null cones, we found that the spherically symmetric 
background spacetime is Minkowskian. We then derived 
a differential equation for a two-tensor, ipAB, whose so- 
lution allows one to determine in a hierarchal manner 
the linear perturbations y AB , C/^, and $(i)- Utilizing a 
representation of the unit sphere with a complex dyad, it 
was shown that the differential equation for i/jab is a wave 



equation for a spin-2 field ?/>.The field ip is corresponds 
to 2r^ , where ^>o is the Newman-Penrose Weyl scalar 
[l7j . We reformulated the hierarchal equations for y AB , 
U^y and $(!) as differential equations for spin weighted 

variables J , U and $, respectively. Since the function 
ip determines all linearized perturbations in vacuo, we 
refer to it as a master function. Under the assumption 
of the existence of a power series of tp in terms of the 
areal distance r at r — 0, we solved the equation for ip 
locally, and subsequently, those for J, U and 3>(i), at 
the vertices. This provided us (cq. ([52^) ) with the lin- 
earized boundary conditions for the Bondi-Sachs metric 
functions in a spin-representation at the vertex in vac- 
uum spacetimes. It also generalizes previously presented 
axially symmetric boundary conditions [22| to the three- 
dimensional case with no symmetries. These boundary 
conditions may be used in numerical simulations to calcu- 
late vacuum space times in the Bondi-Sachs framework 
when the vertex of the null cones is part of the numerical 
grid. 

We employed the boundary conditions for ip in solv- 
ing the wave equation for ip globally by two different 
approaches. In addition, we required the solution for ip 
to be finite at large distances to assure an asymptotically 
flat solution for the perturbations J, U, and In both 
approaches, we represented ip as power series in terms of 
spin-2 spherical harmonics with coefficients depending on 
the retarded time u and the radius r. This allowed us to 
deduce a partial differential equation in terms of u and r 
for the coefficients of this series. 

In the first approach, we imposed a standard ansatz 
of separation of variables to solve the differential equa- 
tion for the coefficients of the series. In the procedure we 
determined an ordinary inhomogeneous differential equa- 
tion (eq. (|59"|0 that is most generally solved by a finite 
spectral series using polynomial coefficients with mod- 
ified spherical Bessel functions of the first and second 
kind as base functions. As the modified spherical Bessel 
functions of second kind are singular at the origin, they 
must be discarded by the regularity conditions at the 
origin, whereas a solution purely depending on the mod- 
ified spherical Bessel functions of first kind obeys this 
regularity condition. Although, this solution (eq. 
for the coefficients of the spectral series of ip is regular at 
the origin, we discarded it, because it diverges exponen- 
tially as r tends to infinity, and it would have generated 
a solution for the perturbations that is not asymptot- 
ically flat. Hence, regularity of the Bondi-Sachs met- 
ric at the vertex is not a sufficient requirement to ob- 
tain an asymptotically flat solution of the Einstein equa- 
tion in the Bondi-Sachs framework. Our calculation also 
demonstrated that using a standard separation of vari- 
ables, where the function is decomposed into a product 
of four factors of which each depends on one of the coor- 
dinate x a only, is unsuited to solve the wave equation in 
Bondi-Sachs coordinates, if one requires the solution to 
be regular at the vertex and asymptotically finite. 

In the second approach, we transformed the second- 
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order wave equation for ip into a inhomogcncous first- 
order transport equation (cq. (|66p ). The inhomogenity 
of this equation vanishes for the lowest (I = 2) spin- 
2 harmonic. The corresponding homogeneous transport 
equation has the characteristic surface u + 2r = const. 
Since the wave equation for ip and the transport equa- 
tion are related by a linear integral transformation, the 
solution of the wave equation for ip and the correspond- 
ing transport equation have the same characteristic for 
the lowest multipolc. Using a polynomial ansatz (eq. 
(|67[l) that incorporates the characteristic information of 
the wave equation, allowed us to find a solution for all 
multipoles for the master function and linearized per- 
turbations (eq. ([71]) ) that are regular at the origin and 
asymptotically flat. This solution is referred to as SPIN- 
2, because it represents spin-2 waves propagating on a 
Minkowski background spacctimc. SPIN-2 has some ad- 
vantages in regard to other linearized solutions in the 
Bondi-Sachs framework found in the literature. First, 
SPIN-2 is given by simple rational expressions in terms of 
the spin-weighted quantities of the Bondi-Sachs metric. 
Second, it describes all multipole of the perturbations, 
while other solutions give only their lowest multipoles 
[HI, Ojl or require an elaborated procedure to generate 
those multipoles fl2l.[bl|. 

For the SPIN-2 solution, we calculated the Weyl scalar 
^4 at null infinity using the formalism of j24j. For ped- 
agogical reasons, wc also summarized the formalism of 
[24| by pointing out the four most important steps in 
obtaining a simpler formula for "J 4 at null infinity in lin- 
earized gravity. This simple analytical expression (eq. 
(I116P 1 ) for ^4 at null infinity and the explicit form of the 
perturbations J, U, and (cq. CLE}) make SPIN-2 
an ideally suited testbed solution for simulations in the 
Bondi Sachs framework and to test numerical wave ex- 
traction methods at null infinity, e.g. (20j describes the 
most recent process of such simulations (for others see 

!)• 

In the future, we plan to test the stability of the SPIN- 
2 solution against small perturbations, to investigate its 
physical reliability. It also would be interesting to utilize 
SPIN-2 to study quadratic perturbations with respect to 
a Minkowskian background spacetime. 



Appendix A: Ricci tensor contributions for the main 
equations 

1. Contributions at O(e ) 

The non- vanishing Ricci tensor contributions of the hy- 
persurface and evolution equations are at O(e ) 



/Vrr , 

r ar 

„ ?r d r 



(Al) 

-@(q AB ) 

(A2) 



where SH^ab) — 2 is the Ricci curvature scalar of the 
unit sphere. 



2. Contributions at 0(e ) 

The relevant Ricci tensor contributions for the hypersur- 
face equations at 0(e) are 



TZr 



(A3) 



2r z 



r*e- 2 e° qAE U E 1)jr 



3^(1 



n EF T> ^ 

"2 9 U El A F,r 



(A4) 



1 



and those for the evolution equations 



-2e 2f3o D A D B (3 {1) + qAEDB^U^ 



-qABD E 



r 2 JJ E 



(A6) 



VIII. ACKNOWLEDGMENT 

It is a pleasure to thank E. Miiller, J. Winicour, A. 
J. Penner, L. Lehner, P. Jofre Pfeil for comments on the 
manuscript and discussions. I also thank B. Carter for 
providing me the proceedings of ref. [U and j34|. Fi- 
nancial support is appreciated from the Max Planck So- 
ciety, from the Collaborative Research Center on Grav- 
itational Wave Astronomy of the Deutsche Forschungs- 
gemcinschaft (DFG SFB/Transregio 7), the Observatory 
of Paris, and the CNRS. 



Appendix B: The master equation and the fiat-space 
scalar wave equation 

In this appendix, it is shown how the master equation 
relates to a flat space wave equation of scalar field with 
spin weight zero. This offers a comparison to the ap- 
proach of [l4| , where the perturbations are generated by 
a spin-0 fields. 

The homogeneous flat space wave equation for a spin-0 
field h is 

= Dh, (Bl) 



16 



where Oh :— ri ah V a Vbh is the d'Alcmbcrt operator, 
which reads in outgoing polar null coordinates 

r 2 Uf = 2r(rh) >ur - (r 2 h ir ), r - cM. (B2) 

If we commute the eth and eth bar operator in the 
master equation ([30a|, we obtain 



= -r 



r 4 (2V>,„ - W) 



3V 



(B3) 



Setting tjj := d 2 F, where F has the spin weight zero, and 
inserting this new definition into (|B3p yields 



r 4 (2F u -F r )l -(33-2)F 



(B4) 



where we again commuted the eth and ethbar operators 
to factor out 3 2 . Eq. (|B5|) implies for F the following 
differential equation 



r 4 (2F„ - f» - (33 + 2) F 



(B5) 



To find how (|B5j) relates to ()Blj) we introduce a further 
spin-0 field F = r n f and insert this definition into (|B5|) 
which yields after dividing out r n 

= 2r(rl r ), u -(r 2 f, r ), r -mf + 2r(n + 3)f, u 

-2r(n+l)/, r -(n + 2)(n + l)/. (B6) 



with an additional damping term — 4r/ iU . Equation (|BT|) 
has also been obtained by [13 in a different approach. 
For n = — 2, eq. (|B6[) becomes 



2 Df 



-2r(f,u + f,r) 



(B8) 



which is a wave equation with a damping and a force 
term. This equation has the disadvantage to (|B7j) that 
it does not reduce to the Laplace equation for time- 
independent fields. Therefore (|B7[) is preferable to (|B8[) . 
and we conclude that the master function ip is related to 
the spin-0 field / via ([B7]) and 



1 



3 2 / 



(B9) 



1> 

Based on (|B7[) and (|B9[) . we now sketch an alternative 
spin-0 approach to the one given in [l4j . Let a and Z be 
spin-0 field that are related to the perturbations J and 
Z via 



J := 
U := 



2q 



3Z 



(BlOa) 
(BlOb) 



Inserting (|B10|) into ([30b ]) -([30d ]) while using (|B9[) allows 
us to deduce three equations 



It can be seen that the first three terms in (|B6j) cor- 
respond to r 2 Df which is the principle part of the flat 
space scalar wave equation in outgoing Bondi-Sachs co- 
ordinates. The other additional terms are indicate that 
(|B6jl is a a quasisphcrical wave equation. In fact, it is 
not possible to set all these terms to zero for any num- 
ber of n. The linear part in (|B6[) vanishes if n = — 1 or 
n = —2, i.e. in this case there are no restoring forces. In 
particular, if n = — 1 and if / is time- independent then 
(|B6[) is the standard Laplace equation in spherical coor- 
dinates. Whereas if / is time-dependent, we have the 
wave equation 



= -4rf. 



(B7) 



( ra ), 



f 



r A Z r 



(33 + 2)f 



$ r 



(33 + 2)a + — (7 



(Blla) 
(Bllb) 

,(Bllc) 



which can be used to the determine the perturbations 
J, U, and $ in the following procedure: (1) solve the 
damped wave equation (|B7[) ; (2) integrate cqs. (IB1 1|) 
according to the given hierarchy to find <I>; and (3) use 
(|BlO|) to obtain J and U. 
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